Imaging of critical correlations in optical lattices and atomic traps 
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We consider critical space-time correlations in trapped utracold atoms and propose a method for 
their detection by the interference with a reference Bose-Einstein condensate. An important infor- 
mation about universal properties can be obtained despite non-uniformity and small sizes of atomic 
clouds. The advantage of this method is that the correlation properties of the interference pattern 
are robust with respect to expansion and are not, practically, affected by microscopic structure of 
underlying lattice including structural disorder. The non-destructive scheme allows measuring the 
space-time correlators of arbitrary order. 
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I. INTRODUCTION 

Experimental realization of the Mott insulator (MI) 
superfluid (SF) transition in optical lattices |jj opened 
new directions in studying strongly correlated states of 
matter. Mott insulator of single-component bosons rep- 
resents a simple and yet very rich strongly interacting 
system. Several other strongly correlated phases have 
also been predicted recently [M 0, 0> HI • The impact of 
disorder on the transition is another area of great inter- 
est 0, which can also be addressed in optical lattices 
(OL) 0. Despite these quite exciting developments, the 
fundamental aspect of the many-body physics — critical 
behavior upon approaching a continuous transition point 
— evades experimental determination in optical lattices. 

The critical behavior is characterized by the divergence 
of the correlation length £ as the system approaches the 
critical point from either side of the phase transition. On 
the MI side, £ can also be viewed as the coherence length 
controlling exponential decay of the phase correlation. 
On the SF side, the system acquires a global average 
phase so that the coherence length is determined by the 
system size L. In this case, the correlation length £ gives 
the length scale of fluctuations on top of the established 
off-diagonal long range order. Thus, £ is finite and coin- 
cides with the healing length (the length scale on which 
the order parameter responds to a local perturbation). 
This divergence of £ is the key feature of the hyperscal- 
ing paradigm of the continuous phase transitions. The 
exponent v characterizes the singular growth of £ as the 
system is tuned into the critical point and uniquely rep- 
resents the universality class of the transition. 

Coherence properties of Bose-Einstein condensate 
(BEC) have been demonstrated in the seminal experi- 
ments [H 0, ^j. In contrast, no systematic studies of 
the critical properties were reported so far. Relatively 
small and non-uniform samples as well as the specificity 
of the imaging methods do not allow direct measurement 
of divergent £. As discussed in ref.pjj, a presence of 
the parabolic trapping potential requires fulfilling spe- 
cial conditions in order to observe critical behavior. For 
some quantities even critical exponents can be modified 
because of the trapping potential. In recent Monte Carlo 



studies 0, it has been shown that practical real- 
ization of the critical behavior requires suppression of 
the parabolic potential and using, instead, quartic traps 
so that the central region of highest density can well be 
approximated by uniform box which is the standard as- 
sumption of the finite size scaling approach (FSS) (see in 
Ref.pj) developed for macroscopic and uniform (apart 
from sharp boundaries) samples traditionally studied so 
far. Thus, consistent application of the FSS paradigm in 
the cases of trapping potential requires separate study. 
In the present paper, we will focus our attention on the 
problem of measuring the correlators of the atomic field 
in OL. Clearly, a solution of this problem is a necessary 
prerequisite for a consistent study of the criticality in OL 
and traps. For sake of concreteness, we shal l ign ore the 
issue of the shape of the trapping potential 0, 0, 0] 
and follow the standard FSS approach [T3 |. 

The absorptive imaging techniques used in. e.g., 
Ref.|l5j| in their present form do not leave much room 
for improvement. Indeed, the far zone imaging of the 
cloud released from an optical lattice for the observa- 
tion of MI-SF transition relies on analyzing the Bragg 
interference pattern. The information on the superfluid 
off-diagonal correlation length is encoded in the width of 
the Bragg peaks 0] . This information can easily be dis- 
torted due to many external causes. Furthermore, this 
imaging method becomes virtually useless for the study 
of superfluidity in the presence of purposely imposed dis- 
order [2|, because rather the lattice disorder than the 
subtle many-body correlations will immediately broaden 
the Bragg peaks. Thus, it is very desirable to have a 
method which will allow a real space imaging of the long- 
range off-diagonal correlations regardless of the underly- 
ing structure of the optical lattice. 

We propose such a method of reference interference 
imaging in which the tested cloud is allowed to inter- 
fere with some reference BEC. Here we consider the sim- 
plest case of one-component bosons undergoing MI-SF 
transition. The question to what extent the proposed 
method can be applied to more complex situations of 
several species is left out for further study. Our idea 
stems from the observation that quasi-coherence in the 
MI phase exists on sizes R such that a; <C R < £, where 
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<x\ stands for the lattice constant. This view is closely 
related to the idea of quasi-condensate introduced, first, 
by Popov 0| and further developed in refs.0,0] with 
respect to weakly interacting Bose gas undergoing a tran- 
sition from a normal state to the BEC. These two seem- 
ingly very different systems can be mapped on to each 
other in the critical region with the help of the quantum- 
to-classical mapping [20j so that the critical behavior of 
the quantum MI-SF transition (at temperature T = 0) 
in d-spatial dimensions is equivalent to that of the Bosc- 
Einstein condensation in D=d+1 dimensions (at finite 
T). In other words, both systems belong to the same 
universality class U(l) characterized by identical criti- 
cal exponents. Thus, studying the critical behavior can 
equally be done either in optical lattices at T = in 
d = 2 geometry or in the d = 3 traps close to the BEC- 
formation temperature |2l| . Fin ding these exponents is 
a very challenging experimental 22] and computational 
task |23. 

It is also worth noting that both transitions are char- 
acterized by a large critical region, which makes it quite 
promising to study the criticality in the trapped atomic 
systems. In some systems such as traditional supercon- 
ductors, mean field behavior is predominant everywhere 
except for a very small region (as small as 10 -15 of the 
transition temperature) about the critical point where 
true critical behavior prevails. For a weakly interacting 
boson gas in free space, the relative width of the critical 
temperature interval around the condensation tempera- 
ture T c is given by the mean density n and the scat- 
tering length a as k 10 2 n 1 / 3 a [ljj. Thus, even for the 
gas parameter as small as na 3 = 10~ 6 , the critical re- 
gion is comparable to T c . Similarly, in the case of the 
zero temperature MI-SF transition in d=2 dimensional 
optical lattices, the critical region width of the quantum 
transition is given by Su/u c rj 1/n, where h is the mean 
occupation number per lattice site, and u = U / J with U 
and J standing for the onsite interaction and the tunnel- 
ing amplitude, respectively |24|. For n = 1, the relative 
size of the critical region is unity. We therefore expect 
that the critical relation £ ss ai/\8u/u c \ v may be used to 
estimate the coherent length reliably down to very small 
length scales. For example, in order to observe £ f=a 10a/, 
one needs 8u » 0.03u c = 0.7 in d=2 square lattice. 

When the correlation length becomes mesoscopic, a 
particular region of linear dimension ~ £ can be viewed 
as a domain of coherent SF. Such a finite size SF is char- 
acterized by uniformity of the classical phase. If the op- 
tical lattice is removed and the cloud is allowed to over- 
lap with some reference BEC cloud moving coherently 
at momentum k r S> l/£, an interference pattern will be 
formed. This pattern will consist of domains of typical 
sizes ~ £, with each domain featuring a typical interfer- 
ence pattern of parallel lines separated by the distance 
2n/k r 0. Pronounced phase shift (or distortion) in the 
parallel lines will occur over the length scale of £. 

This proposal has a very close analogy in quantum 
optics j25(: Measurements of the phase of an optical 



beam are generally performed by sending the beam onto 
one input arm of a beam-splitter, and the (strong) ref- 
erence beam onto the other input arm. The intensity 
and noise in the two output arms carry information on a 
given quadrature of the optical beam, and on its variance. 
In particular, squeezing experiments consist in observing 
a reduced noise on a quadrature. In the matter waves 
field the interference between two condensates has been 
employed for identifying vortices |2(|. 

It is important to note that the relative phase between 
the optical lattice cloud and the reference BEC varies 
from shot to shot in an uncorrelated way, so that no in- 
terference pattern will appear when averaged over many 
realizations. However, as we will show below, the combi- 
nation of the interference with the measurements of the 
noise correlator of the Hanbury Brown & Twiss (HBT) 
type [13, |2iJ |3(j , which gives the density-density cor- 
relator at different spatial locations, can provide the in- 
formation on diverging correlations. Related methods 
have been recently proposed to determine long-range or- 
der in fermionic systems [3lj . Very recently, the noise 
analysis of the interference between two low-dimensional 
condensates has been proposed for detecting the specific 
fluctuating behaviors [33 • We will also discuss that the 
reference interference imaging can be adopted for deriv- 
ing a full space-time correlator, which carries the most 
complete information about a system, including excita- 
tion spectrum and the dynamical scaling. Furthermore, 
the non-destructive variation of the method employing 
outcoupling atomic beams allows, in principle, measur- 
ing higher order correlators - n-body Green's functions. 



II. CORRELATION PROPERTIES OF THE 
INTERFERENCE PATTERN 

In the following subsections, we, first, discuss the 
Bragg interference of the lattice atoms alone in order to 
emphasize the limitations of the Bragg method. Sec- 
ond, we will consider the interference between the lattice 
atoms with a reference BEC so that the advantage of 
using the reference interference imaging will become ev- 
ident. The analysis of a simplified scheme, in which the 
reference cloud is taken as uniform, will be followed by a 
discussion of a more practical setup later in Sec lIVI 



A. Density imaging of expanding cloud released 
from a lattice 

Earlier experiments focused on expanding density pro- 
file of atoms released from a trap or a lattice, and here we 
give a brief discussion how this may or may not reflect the 
superfluid correlations in the lattice. In general, density 
profile n(x, t) of a cloud undergoing free expansion is de- 
termined by distribution of momenta carried by particles. 
In other words, n(x, t) represents spatial Fourier trans- 
form of the one-particle density matrix with momentum 



3 



q = x/t. (Here and below, unless stated otherwise, we 
employ atomic units in which atomic mass m and the 
Planck constant h are unity). 

In the absence of the lattice potential, when £ is the 
shortest length scale which defines the largest momen- 
tum l/£, the expanding cloud size R at long times would 
be uniquely determined by £ as R — t/£. In the presence 
of the lattice, higher spatial harmonics impose a back- 
ground which has nothing to do with the critical behavior 
and which complicates the analysis. 

Onsite packets of Wannier functions Wi(x — Xj), where 
x* is the position of site i, contribute the fastest compo- 
nents into the expansion. Their expansion speed is deter- 
mined by the uncertainty principle V w ~ l/a w , where a w 
is the width of the initial Wannier function. Thus, a typ- 
ical size of the packet becomes R w « t/a w after a brief 
initial transient on the order of t w w o? w . Expansion of 
each Wannier packet introduces an additional phase fac- 
tor exp(iA) given by the classical action of a free particle 
A = (x — x') 2 /(2t), where x' and x are the initial and 
final positions respectively. Therefore, the result of free 
expansion of a Wannier function can be represented as 

Wi(x - x 4 ) = j^e 1 ^ 1 J dx'e^^^^Woix'). (1) 

Here, Wo denotes the Wannier function at the initial 
time, and we have neglected a phase of x' 2 /2t which is 
much smaller than unity after the transient when x' is 
taken within the range a w of the initial Wannier func- 
tion. The factor of t~ 3 / 2 takes care of the decrease in the 
magnitude as a result of the expansion. 

Simplicity is achieved after a relatively long expan- 
sion time, t/a w >> L, where L is the linear scale of 
the lattice. By now, the Wannier function has expanded 
much beyond the original size of the lattice. In this 
case, the phase term Xjx'/t is small, and one arrives at 
the expression traditionally used for the time of flight 
imaging^ W t (x - x 2 ) = ^ 3 / 2 exp(i(x - x;) 2 /(2t))Wo(q), 
where H^o(q) stands for the Fourier transform of Wo(x) 
at q = x/t. The expression for the evolution of the 
Wannier function can be employed in the field operator 
for the lattice atoms 

^(x,t) = ^W t (x-x i )6 i , (2) 

i 

where 6, is the destruction operator for the Wannier state 
on the ith site of the lattice. Apart from an unim- 
portant numerical factor, the imaged density n(x, t) = 
(?/>/(x, t)V>z(x, t)) becomes 

n(x,t) = I|^ (q)| 2 5(q), (3) 

x 2_ x 2 

5(q) = ^e^-^e'V^}. (4) 

ij 

Further simplification can be achieved in the far 
zone, when t » L£. The phase factor e % at can 



be then ignored. It allows to relate n(r) directly to 
the Fourier transform of the lattice correlator S(q) = 
2y exp(iq(rj — v.j)){b\bj). This expression indicates 
that, as long as the OL is regular enough, the Bragg 
structure of peaks will be seen, with their width deter- 
mined by 1 /£ (or eventually limited by finite size effects in 
the critical regime where £ > L). However, if the lattice 
is not regular, the correlations in {b\bj) will be strongly 
obscured by the factor exp(iq(r^ — r,-)). In what fol- 
lows, we will show that this limitation can be overcome 
by measuring the interference pattern with a reference 
BEC. 



B. Interference with a reference BEC 

For illustration, we imagine that immediately after the 
lattice was instantaneously removed, a coherent matter 
plane wave of the form exp(ik r x)/v / f2, normalized by 
some large volume f2, is sent onto the lattice and is 
made superimpose with the lattice atoms (more details 
are given in Sec. III). The result of a single-shot absorp- 
tive imaging is given by the state average of the total 
density n(x) = (S\ip^ (x)ip(x)\S) , with the field operator 
represented as 

^ (x) = w *(* - *i) b i + 6 -/' ( 5 ) 

i v 

where b r f is the destruction operator in the reference 
BEC state. The total density contains an interference 
term proportional to 

(S\bl f bi\S) ttT/N^he^, (6) 

where N r f stands for the total number of atoms in the 
reference cloud, whereas h and characterize the coher- 
ent nature of the lattice atoms in the state \S). 

According to the theory of critical phenomena, the de- 
struction of the SF phase occurs because of large fluc- 
tuations in the phase ipi while the amplitude n remains 
relatively constant. For this reason, one can talk about a 
quasi-condensate for the lattice atoms through the crit- 
ical region of the phase transition. One may view h as 
the quasi-condensate density per site, and ipi as a classi- 
cal phase with the same significance as that of a classical 
wave. Indeed, for £ d h ^> 1, the state \S) can be rep- 
resented in the coherent state picture specified by these 
parameters |33). 

Although these phases vary randomly from shot to 
shot, they may be regarded as fixed in a single shot inter- 
ference measurement 8] . Spatial dependence of the clas- 
sical phase contains crucial information about the state 
of the system. Deep in the SF phase, <fi is uniform in 
space. Upon approaching the critical point, it fluctuates 
but retains a uniform part, which becomes zero at the 
critical point. On the MI side of the phase transition, 
the uniform part of the phase does not exist anymore in 
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the sense that the spatial average of exp(i<£>i) vanishes in 
thermodynamical limit. The length scale of phase fluctu- 
ations is the correlation length £, which diverges upon ap- 
proaching the critical point from either side of the phase 
transition. 

How can one observe the spatial behavior of the clas- 
sical phase? The interference fringes are given by the 
cross-term in the density 

Sn(x) = A(x) e ~ ikrX + c.c, 
A(x) = ^Jfi W (x - XJ )e^ , (7) 

3 

where n r f = N r f/Q is the reference cloud density [34| . 
Assuming that the correlation length is much larger than 
the lattice constant, we can replace <fij by y>( x ) i n the ex- 
pression of the amplitude EqJJJ), where <p(x) is a smooth 
interpolation of the classical phases on the lattice sites. 
The interference fringes are then described by a sinu- 
soidal function of k r x — ip(x) with a modulating ampli- 
tude given by the magnitude of \fn r ffi^- Wq(x — xj). 
Such amplitude exhibits strong spatial oscillations on the 
length scale of the lattice constant. Already at expansion 
times £ 2 ^> t ^ La w , the fast oscillations have almost 
completely been smoothened out (what remains of them, 
can still be filtered out by means of Fourier analysis). At 
the same time the phase domains of the size ~ £ are still 
well separated from each other and feature the interfer- 
ence pattern of parallel stripes separated by a distance 
27r/fc r . 

In practice, the role of larger expansion of the lattice 
atoms after the release of the optical lattice must be con- 
sidered, because it is difficult to carry out the measure- 
ment in the short time before expansion. However, the 
expansion does not change the correlation length. Elim- 
ination of the OL transforms the cloud into a weakly in- 
teracting system so that the subsequent expansion can be 
considered essentially free. Indeed, the correlation prop- 
erties are determined by a typical width of the momen- 
tum distribution in the system. Release of the OL essen- 
tially freezes in this distribution, which implies "conser- 
vation" of £. Therefore, the qualitative features discussed 
above still apply to the interference pattern between the 
expanded lattice- atom cloud and the reference BEC. 

After smoothening out the short range modulations, 
the interference pattern consists of parallel stripes sep- 
arated by a distance determined by the wavevector k r . 
Spatial variation of the classical phase gives rise to distor- 
tion of the stripe positions on the scale of the correlation 
length. The stripe pattern is visible in a single-shot mea- 
surement provided each domain contains a large number 
n£ d 3> 1 of atoms so as to overcome shot noise. Statisti- 
cal analysis on the distortion of the interference pattern 
can then yield crucial information on the spatial corre- 
lation of the classical phases: CV, = (e lVi e~ ttpj ), where 
(. . .) stands for average over the whole lattice for a fixed 
distance between the sites. Such a correlator is insensi- 
tive to the random variation of the phases from shot to 



shot, and can thus be enhanced by a further average over 
many runs of the experiment. What do we expect on the 
quantitative behavior of this correlator? 

On the MI side of the phase transition, the correlator 
decays to zero as Cy ~ exp(— |x.;— Xj\/^) for |x^— Xj| > £. 
For |xj — Xj| < £, the correlator exhibits a power law 
decay typical for the critical point ~ l/|xj — Xj\ 1+V with 
some universal exponent rj w 0.038 in 2d lattice across the 
MI-SF transition and in a 3d gas undergoing the (finite 
temperature) BEC transition. 

On the SF side of the phase transition, the correlator 
retains a finite value Coo characterizing the off diagonal 
long range order as |xj — Xj\ — > oo. This order param- 
eter vanishes as a power law CLo ~ \u — u c \ 2 ^ , with an- 
other universal exponent j3 |14j, as the critical point is 
approached. The decay from Cu = 1 to Coo << 1 is 
described by the correlation length £ on the SF side. 

The critical point is practically reached when the cor- 
relation length exceeds the system size L. In this case, 
the order parameter behaves as Cy, ~ L~ 2 PI 1 ' in accor- 
dance with the finite size scaling 14]. This corresponds 
to a power-law decay of the fringes with the system size 
L, which highly contrasts with the standard interference 
situation of two coherent BEC clouds deep in the SF 
phase Q, where no L-dependence is to be anticipated. 

C. Density-density correlations in the interference 
pattern 

At long expansion times t > £ 2 , the phase domains 
start overlapping so that visibility of the fringes with the 
reference BEC is strongly reduced. However, effects of 
the coherence of the initial many-body state can still 
be observed in the density correlations, which can be 
measured by taking the average over many realizations 
|29|. Such experiments are matter- wave analog of the 
Hanbury-Brown and Twiss (HBT) effect j2?J widely em- 
ployed in quantum optics, and important results in this 
direction have recently appeared [2^, [3(j • 

The density-density correlation function is defined by 

if(xi,x 2 ;t) = (n(x 1 ,t)n(x 2 ,t)), (8) 

where the density operator is given by n(x, f) = 
^t(x,f)^(x,f) with 

#M)=^(x,i)+Vv/(x,t). (9) 

The subscripts I and rf refer to the states of the lattice- 
released and the reference clouds, respectively. We as- 
sume as before 

e ik r x 

ip rf (x,t) = -j=-b rf (10) 

(a more practical setup will be discussed in section III). 
For the cloud released from the lattice, the operator 
tpi (x, t) is given by eq. @ ■ 
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The density operator can be split into three terms 
due to: lattice-released cloud n t — ipjipi, the reference 
cloud n r f — ipljiprf, and the interference between the 

two rimt = iplfipi + tpjiprf, respectively. Then, among 
the many terms which result, we are mostly interested 
in the correlation function for the interference term rii n t 
alone. In the center-of-mass r = (xi +X2)/2 and relative 
p = X! — x 2 coordinates, it reads: 

K int (r, p; t) = (rcint (r + p/2) n int (r - p/2)) = 
n rf j^-P {i>]{r + p/2, t) Mr ~ p/2, t)) + c.c. (11) 

This term is directly proportional to the one-particle 
density matrix in OL modified by free expansion. Be- 
low we will show that the integration over the center of 
mass coordinate r results in the quantity Ki nt {p,t) — 
Jdr Ki n t(r, p,t) which represents spatial correlations av- 
eraged over the lattice and which turns out to be invari- 
ant in time and therefore insensitive to the expansion 
|35j ] . In fact, Kint represents the mean of the single par- 
ticle shift operator (exp(ippi)), where pi stands for the 
momentum operator for one particle. During free ex- 
pansion, the momentum of each particle is conserved, 
implying that the mean shift is a constant in time. 

We comment on the strength of the remaining terms. 
The cross correlation between m n t and n; + n r vanishes 
upon ensemble average. The remaining terms are from 
the correlator of ni + n r f. The corresponding terms pro- 
duce either a uniform background or have a different 
spatial periodicity than k r , and can thus be filtered out 
from an actual measurement by means of spatial Fourier 
analysis. The correlation function for the interference 
density can also be identified based on a strength analy- 
sis, if we have a relatively strong BEC for the reference 
cloud and a relatively low density for the lattice-released 
cloud. (The latter scales as 1/t 3 after expansion, which 
can be made weaker by waiting longer.) In this case, it 
is even more convenient to consider the noise correlator 
K'(r,p,t) = (An(r + /9 /2,i)An(r- p/2,t)}, which 
has no cross terms with respect to the three parts of the 
density operator. The fluctuations An = n — (n) will be 
suppressed in the reference cloud as long as the BEC is 
well formed. The noise correlator for the lattice-released 
cloud is small (scales as t~ 6 ) due to the assumed low 
density. Thus, the term l(TT)l. which scales as ~ n r ft~ 3 
turns out to be dominant. In our subsequent discus- 
sions we shall not explicitly distinguish K'(r,p,t) and 

K in t(T,p,t). 

We can thus focus on the space-averaged correlation 
function at time t = 0, which can be evaluated as 

K int (p;t = 0) = n rf e*"-' £ (b\bj) e^ {p ~ Pii? + c.c. 

ij 

(12) 

Here, we have assumed for simplicity a Gaussian form for 
the initial Wannier function: 

M/ (x) = a - 3 / 2 e-* 2 / 2 < (13) 



The quantity (b\bj) may be regarded as equal to the 
quasi-condensate density n times the phase correlator Cij 
discussed in section IB. Let C(p) be an envelop function 
that equals the phase correlator at p = pij and smoothly 
interpolates between the lattice vectors. We expect this 
envelop to vary on the scale of the correlation length £, 
which is taken to be much bigger than the microscopic 
size of the Wannier function a w . Then we can write 

K int (p; t = 0) = n rf e^'P C(p) £ ne~^ (p ~ Psj)2 + c.c. 

ij 

(14) 

The phase correlator envelop of the lattice thus enters di- 
rectly as a factor in the density-density correlation func- 
tion. The other important factor is the plane-wave mod- 
ulation from the reference BEC. The factor from the lat- 
tice sum varies over the scale of the lattice spacing and 
can be easily separated out. In the following analysis, we 
assume that a smoothening procedure has been applied 
to this factor, so that it is replaced by a number pro- 
portional to the total number N of the particles in the 
lattice. 

Having considered the behavior of the phase correlator 
in Sec III Bl we can now summarize our results for various 
regimes for the quantum MI-SF transition in d=2 OL and 
in d=3 traps undergoing finite temperature BEC transi- 
tion. As mentioned above, both systems are equivalent 
in terms of the universality of the transitions. Here we 
assume that the condition of weak gradients is sat- 
isfied in both cases. This condition can be formulated in 
terms of a single length ~ L determining both the gra- 
dients and the system size, so that as L is taken larger 
and larger the system becomes progressively closer to the 
transition and the typical scaling (with L) behavior can 
be observed. In both cases for £ < L, in the phase were 
no symmetry is broken 

K int (p, t) ~ n rf N + c.c. (15) 

In what follows, we will call this phase normal with- 
out specifying if it refers to the MI or to the ther- 
mal gas above the BEC transition. When £ becomes 
larger than L, a system is in the critical regime, where 
K int (p,t) ~ e^P/p 1 ^. The exponent r] 14] is related 
to the exponents (3 and v by means of 2(3 /v = 1 + n. 
This implies, in particular, that measuring the magni- 
tude of the correlator at p w L, which goes like L~ 2/3 / u , 
in progression of L can provide valuable information on 
the exponent f3. These expressions are independent on 
particular details of a system undergoing the phase tran- 
sition and are not distorted by the free expansion. 

In actual experiments, it is the columnar density that 
is usually measured in absorptive imaging. In the case 
of a 3d BEC transition in a trap, the functional depen- 
dence in (|15|) is not changed much. If the columnar 
density is collected along z— axis which is normal to k r , 
the 3d function (|15fl must be replaced by the integral 
n rf Ne lW ' r J dze-^ r ' 2+z ' 2 lt/(r 2 + z 2 )^^ 2 , with r being 
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2d coordinates in the xy plane. As a function of r/£, this 
integral exhibits exponential behavior similar to ea. (|15fl . 
In the critical situation, with the correlator measured on 
distances ~ L, the columnar density interference pattern 
will acquire one extra power of L as ~ l}- 2 P/ v = L~ n . 

The 2d MI-SF transition can naturally be realized in 
a stack of identical and independent 2d lattices, and the 
columnar density can be collected along the z-axis per- 
pendicular to the layers, with k r being along the layers. 
Differently from the 3d case, this simply leads to a mul- 
tiplicative prefactor in (15), so that the correlator scales 
as ~ L~ 2I3 / V = L~( 1+r >), 

Useful results can also be obtained from the Fourier 
transform Q(£,L) = J dpKi n t(p, t) cos(k r p). In 2d OL, 
it must scale as 

Normal: Q(£, L) ~ iV£ d_1 = L d £ d_1 , (16) 
Critical : Q{^,L) ~ NL d ~ 1 ~'' 1 = L 2d ~ 1- '', (17) 
BEC : Q(L) ~ TV 2 = L 2d , (18) 

where the last line is given for the BEC case in which 
the phase correlator in i|15|) must be replaced by just 
a constant. The dimension d should be understood as 
follows: d = 2 in the MI-SF case (which can be realized 
as a stack of independent 2d-layers) and d — 3 in the case 
of the BEC transition in a 3d trap. As we will discuss 
later, these properties are quite robust with respect to 
the clouds non-uniformities (provided these are weak in 
a sense discussed above). Thus, observing how the global 
characteristic Q(£, V) scales with L can provide direct 
information on criticality. More detailed discussion of the 
cases will be given later, where the practical issues such 
as finiteness and non-uniformity of the reference cloud 
are considered. 



III. NON-DESTRUCTIVE MEASUREMENTS 
OF THE SPACE-TIME CORRELATOR 

A very important information about the system 
can be obtained if the in-situ space-time correlator 
/C(x2, r; xi, 0) = (ipj(x2, t)^;(xi, 0)) is known, where r 
refers to the time evolution inside the lattice or the trap. 
Specifically, the excitation spectrum is provided by the 



r— dependence. Obviously, the destructive scheme can- 
not yield such correlator - the method described above 
gives the one-particle density matrix, that is, the one- 
time correlator corresponding to r = 0. 

A very flexible method of non-destructive detection 
|36| of the phase correlations has recently been intro- 
duced. It relies on interference of the outcoupled atomic 
beams. Another method of the space selective non- 
destructive imaging has been emp loyed in combination 
with the single atom detection |37j ]. 

Here, we show that the method of the noise spec- 
troscopy applied to the images of interference be- 
tween two outcoupled atomic beams and two refer- 
ence BEC clouds with fixed relative phases can pro- 
vide the complete correlator /C(x 2 , r; xi, 0). This state- 
ment can be illustrated in the frame of two suc- 
cessive measurements of the interference contributions 
Ai(xi) = VL/(xi)^ o (xi,0) + H.c. and A 2 (x 2 ) = 

^2 r /( x 2)V'o(x2, t)+H.c. taken at the time moments t = 

and t — t, where ^j r y(xi) and V4r/( x 2) re f er to the 
operators of the reference clouds with well defined rel- 
ative phases. If the measurement is performed rapidly 
enough, the field operators ?/; o (xi,0) and ^ (x2,r) de- 
scribing the outcoupled atomic beams are replicas (mod- 
ulo multiplicative factors depending on the details of the 
extraction scheme) of the field operators describing the 
lattice cloud i/>j(xi,0) and ^(x 2 ,r), respectively. 

Following the same logic applied to the case of the 
destructive imaging, as long as the correlation length 
is large (includes many particles in the volume ~ £ d ) 
one can interpret the one-shot two interference images in 
terms of the phases of the two reference clouds <fi r f, ^2rf 
and of the outcoupled beams <^(xi,0), </?(x 2 ,r). Thus, 

A fc (x fc ) = y/ntofn{*kM)j (v(Xk - tk) -' Pkrt{ - Xk)) + cc., 
where n^rf and n(x/ c ,t/ c ) stand for the densities of the 
k = 1,2 reference clouds and of the outcoupled beams 
taken at the time moments tk, respectively. 

The single-shot interpretation can be extended into the 
correlation multi-shot analysis. As the phases of the ref- 
erence clouds are well defined with respect to each other 
(this can be achieved by splitting one BEC), they can 
be taken out from the average describing the multi-shot 
correlations: 



(^(xa)Ai(xi)) = ^ lr/ n 2r/ n(x 1 ,0)n(x 2 ,T)(e l ^( X2 )-^( Xl ) + ^ x -°)-^( x - T )]) = ^(x 2 , Xl )/C(x 2 , r; x 1; 0), (19) 

I 



where the multiplicative factor „4(x 2 ,xi) = 
^ fr ^j T ^ e i( v , 2rf ( X2 )- l p lrf ( Xl )) ig due tQ the refer _ 

ence clouds. Expansion can be described in the same 
way as in the destructive case, in particular it is easy 
to see that it does not affect the center of mass average 



K(p,r) = J rfr/C(r + p/2,r;r - p/2,0). Thus, studying 
it as a function of the spatial and time separations 
will allow extracting crucial information on the critical 
space-time behavior [20j . 

The above scheme can be extended to the case of 
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higher order correlator |38j , i.e. a n-body Green's 
function of the form K n (x±, t\) . . . ;X2„, i2«) = 
(#t (xi , h ) . . . *t ( x „ , t n )V (x„+i , t n+1 ) . . . * (x 2 „ , t 2 n)) 
calculated at the points (xi, t\); ...(x2 n , *2n) in space- 
time. As 2n time points are involved, one needs to 
prepare 2n reference clouds with well defined relative 
classical phases tpk,rfi k = l,2...2n, and to outcouple 
2n atomic beams at the corresponding time moments 
ti, ti...ti n . Each atomic beam should be allowed to inter- 
fere with a corresponding reference cloud, so that single- 
shot images give the amplitudes A k (x k ), k — 1, 2, ... , 2n. 
Then, the multi-shot correlation analysis provides the 
correlator K n (xi , t\ ; . . . ; x 2 „ , t 2n ) ■ 



IV. PRACTICAL ISSUES 

In the previous discussion we have adopted a simpli- 
fied approach which ignores non-uniformity of a reference 
cloud as well its expansion. In fact, under these condi- 
tions, it is hard to achieve statistical independence of the 
both clouds. In reality, a reference cloud of some finite 
size R > L should be created some distance x > R 
away from the OL so that there is, practically, no overlap 
or it is exponentially small so that it can be considered 
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FIG. 1: A typical interference density distribution for the 
noise correlator K'(p,t) (scaled in arbitrary units) for a stack 
of 2d OL layers of size 32x32 (unit cells) as given by eas. 1221 
123 IHJ . The correlator (b\bj) is taken in the form 1251261 . 
with £ — > oo. 



This implies that J d d pKi nt (p; t) = for any t be- 
cause it is simply proportional to the quantum mechan- 
ical overlap of the one-particle wavefunctions of the OL 
and the reference cloud, which, initially, was taken as 
zero (or exponentially small). [The form l|15|l does not 
satisfy this condition due to the approximation of infi- 



nite size of the reference cloud.] This can be seen ex- 
plicitly by noting that J d d pK int (p;t) = J dxi J dx 2 = 
(n(xi, t)n(x2, t)} and that the integral of the total density 
/ dxn(x, t) = J dx?p}(x, i)V>z(x, t) + J dxtpl(x,t)ip r (x,t) 
because the interference part J dx.ip](x, t)ip r (x, t)+H.c. is 
essentially zero due to orthogonality (or negligible over- 
lap) of the one-particle states of the tested and refer- 
ence clouds. However, the Fourier harmonic Q(£, L, t) = 
J d d pK int (p; t) cos(k rj o) is not an integral of motion and 
does depend on time, reaching its maximum at the closest 
approach of the centers of mass of the two clouds. 

Here we shall consider a situation when both clouds 
are released the distance x apart at t — with some rel- 
ative momentum k. r toward each other to achieve maxi- 
mum overlap of the centers of mass at the time moment 
t = x a /k r (in atomic units m — 1, ft = 1). During the 
free flight period, the expansion due to quantum pressure 
deforms the clouds and the resulting image. Accordingly, 
the correlator Ki nt (p; t) cannot be considered as a con- 
stant of motion anymore, and the time-dependence must 
be calculated explicitly. Further image analysis should 
be based on the fit procedure discussed below. 

Instead of the simplified form (|10f) . we now take the 
reference cloud wavefunction at t = as 



Vv/(x,< = 0) 



1 



ik r (x— x Q 



brf, (20) 



(V^Ro) 3 / 2 

where x D stands for its initial position. Then, after time 
t > of free expansion it evolves into 



Vv/(x, t) 



Rl 12 



2k r (> 



.(*))- 



.or 



7T 3 / 4 (i?2+^)3/2 



2 («s+«) 5 



rf, 



(21) 

where x Q (t) = x c + k r t. This form must be used in eq.@ 
together with eq.J5J and eq.JTJ in the limit R a ^> t/R a ^> 
a w . Then, integrating over the center of mass coordinate, 
as was done in derivation of ea. 1)15(1 . we obtain the 
correlator 



K int (p;t) 



n rf 
i 3 



^2(blbj)e s »coB((pij), 



where the following notations are introduced 
R 2 



(22) 



(23) 



(fij = k r p + -(R 4J + x Q - tk r )(p - p^), (24) 

and = (x^ +Xj)/2 (the OL spacing we take as a unit 
of distance), with the initial positions of the lattice 
sites counted from the origin located in the middle of the 
lattice. The overall numerical factor in ea. (|22|l is omit- 
ted. For demonstrative purposes we consider a model 
form of the lattice correlator on the insulator side of the 
transition as 



(b\bj) 



Gi 



exp(-Gy) 



|x, 

~2 



2 i 



X 



L 2 



(25) 
(26) 



8 



where the term 



L 2 



reflects the non-uniformity in the 



trap and rj 0.04. In the BEC state: 



exp(- 



L 2 



-)• 



(27) 



The spatial profile of the correlator Ki nt (/?; t) is presented 
on 2d columnar density plot FIG. ^ for the case £ ^S> L in 
cas. l|25l26[) and for the moment of the maximum overlap 
t = x /k r . Typical profiles of the noise correlator for the 
BEC case are shown on FIG. E and FIG. [2 




FIG. 2: (color online) The profile K'(p x = 0,p y = y,t) 
(arbitrary units) for a stack of 2d OL layers in the BEC state 
for two sizes L — 16, 32. The horizontal axis is the OL y- 
coordinate in units of ai. The typical pyramidal shape is due 
to the non-uniformity of the OL cloud. 



The integral Q(£,L) becomes 



Q(f,L) = J d d pK mt (p;t) cos(k r p) 



(28) 



The contribution ~ (R^ + x Q — tk. r )(p — pij) to the 
phase factor ipij in ea. l|22|) can destroy the correlations 
imprinted in (b\bj^. However, if R a ^> L, then, at the 
maximum overlap of the centers of mass of both clouds 
(at t = x /k r ) this term can be estimated as ~ L/R a <C 
1, and, thus, becomes insignificant. In what follows we 
choose the following parameters R — 4L, x Q — 8.5L, t = 
8L. Accordingly, k r = x a /t « 1.1 , so that the wave- 
length 2ir/k r ?s 6 in units of the lattice spacing. 

The quantity Q(£, L) is plotted for the case (|Trj|) in 
FIG. H On FIG. E the quantity Q(£,L)/L d , d = 2 is 
represented for the critical correlators (£ = oo) ) (|25|) for 
two values of the critical exponent r/ — and rj w 0.04 
in order to demonstrate how a small value of r\ can be 
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FIG. 3: (color online) The same system as in FIG. [5] with 
the profile taken along p x — x, p y = 0. 
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FIG. 4: (color online) The integral Q(£,L)/L d , d = 2 as a 
function of the system size L (in units of ai) for three values 
of the correlation length £ = 4, 10, 20 shown on the graph. 
The curves are normalized by Q(£,,L = 4)/4 2 . 



resolved. Direct fits of these dependencies from ea. <|17|) 
give 77 = 0.00 + -0.001 for the case 77 = and rj = 0.036 + 
—0.001 for the case i] = 0.04. For a comparison, the case 
corresponding to the BEC is presented as well. Thus, 
this graph clearly demonstrates the drastic differences 
between critical and fully coherent clouds. Obviously, the 
scaling dependences are not obscured by the clouds non- 
uniformity because a single size L <C £ determines the 
lattice size, the dimension of the reference cloud R ~ 
L and the extent of the non-uniformity as implied by 
eq.lHIj!. 
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tions in optical lattices and traps within the concept of 
the finite size scaling. Its essential component is the 
noise-correlation analysis of the pattern obtained from 
the interference with some reference BEC cloud. The in- 
terference term turns out to be, on one hand, sensitive 
to the diverging critical correlations, and, on the other 
hand, essentially insensitive to the underlying structure 
of a lattice, so that this method can be also used to study 
the critical behavior in lattices with imposed disorder. 
The finite size scaling concept can be effectively applied 
for extracting the critical exponents from the integral 
quantity L). The non-destructive scheme employing 
several reference clouds with well defined relative phases 
is suggested as a tool for obtaining many-body, multi- 
time correlators. 



FIG. 5: (color online) The integral Q(£,L) normalized by 
Q(£, L — 4) as a function of the system size L (in units of a;) 
for three cases: 1) labeled as 77 = 0, represents the correlator 
1251261 with 77 = 0; 2) labeled as rj = 0.04, the same as 
in the case 1) except that r\ — 0.04; 3) labeled as BEC/25, 
corresponds to the BEC state as given by eg. 1271 and reduced 
by the factor of 25 to fit into the plot. 
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